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INTRODUCTION 

A  number  of  transport  phenomena  are  reasonably  well  described  by  a 
combination  of  a  flux  equation  and  the  equation  of  continuity.   Among  these  are 
the  diffusion  of  matter,  the  flow  of  heat,  and  the  flow  of  water  in  porous 
media. 


In  the  flow  of  water  in  porous  media,  the  flux  equation  that  applies 

»  law  {2) A       The  comb 
equation  leads  to  the  equation: 


3/ 

is  Darcy's  law  (2).    The  combination  of  Darcy's  law  and  the  continuity 


If  =  V  .  (K(9)   VH)    t  tl] 

Where  0  is  the  volume  water  content,  H  is  the  hydraulic  head,  K  (0)  is  the 
conductivity  as  a  function  of  water  content,  and  t  is  time  (2).   There  are  a 
number  of  assumptions  inherent  in  the  use  of  equation  1.   It  is  not  the  purpose 
of  this  report  to  examine  these.   Items  (2  and  3)  in  the  Literature  Cited  con- 
tain additional  material  on  this  aspect  of  soil  water  flow. 

Equation  1  can  be  modified  as  follows: 

|.  V  .  (D(9,   V  8)+^M  [2] 

where  the  soil  water  diffusivity,  D  (0),  is  given  by: 

D  (0)  =  K  (0) 
d0 
dh  [3] 

The  pressure  head,  h,  has  been  assumed  to  be  a  single  valued  function  of  0. 
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If  within  the  flow  system,  consumption  or  production  of  water  occurs, 
this  may  be  treated  by  adding  another  terra  in  equations  1  or  2  so  that  they 
become : 

§f  =  V  .  (K(9)  VH)  +  A  (x,  y,  z,  t)  [h] 

=  V  .  (D(0)  V9)  +  ^^1  +  A  (x,  y,  z,  t)  [?] 

where  the  function  A  (x,  y,  z,  t)  will  be  denoted  as  the  "source  term". 

It  is  the  purpose  of  this  report  to  describe  a  power  series  method 
of  solution  of  a  one- dimensional  form  of  equation  2  where: 

(1)  The  flow  is  in  the  horizontal  direction  and 

(2)  The  diffusivity  is  either  an  exponential  or  a  linear  function  of 
water  content. 

Then  equation  2  becomes: 

I  -  i  <D  <*>  t>    .  [6] 

The  Exponential  Diffusivity  Function. 

Experimental  evidence  gives  some  indication  (1)  that  the  diffusivity 
function  for  many  soils  can  be  characterized  by  an  exponential  function  of  the 
form 

D  (0)  =  OLe   /S9,   <*  >  0,  J3>   0  [7] 

where  0C  and  /3   are  parameters.   Using  this  function,  we  shall  solve  the  follow- 
ing boundary  value  problems: 

H  =  i(«^9f)>        °<x<OO)t<0,  [8] 

9   (0,    t)  =  9k       (k  =  1,  2),    t  >  o  [9) 

8    (x,    0)    =   9  x  >  0, 

0  <  92  <  9Q  <  Bj         ,  [10J 

The  problems  corresponding  to  k  =  1  and  k  =  2  will  be  termed  the 
Horizontal  Infiltration  Problem  and  the  Horizontal  Outflow  Problem,  respective- 
ly. 

Despite  the  similarity  of  the  two  problems,  the  corresponding 
solutions  exhibit  distinctly  different  characteristics. 


By  introducing  the  dimensionless  parameters  C  =   x/L,  7"=  C*-t/L  and 
making  the  substitution: 


9   =h  InV,  [U] 


equation  8  becomes: 


>    •  [12] 


L  is  an  arbitrarily  chosen  characteristic  length  different  from  zero. 
Solutions  of  equation  12  of  the  form £7  (£  ,  J*  )  =  f  (79)  will  be  sought,  where 
77  =  C   hf2T  .   Then  equation  12  reduces  to  an  ordinary  nonlinear  differential 
equation  and  the  preceding  problems  become: 

f  (71)  f"  (7J)  =  -  !{  V   (7|)  [13] 

f  (0)  =  e^9*,   (k  =  1,  2)  [Ik] 

f  (7[)   =  e  ^9o  [15] 


0  <  02  <  9     <  0 


[16] 


The  resulting  two-point  boundary  value  problems  pose  difficulties 
because  the  condition  of  equation  15  is  not  easy  to  satisfy.   Fortunately,  the 
differential  equation  13  and  the  initial  conditions  are  such  that  each  of  the 
family  of  curves  satisfying  the  latter  tends  to  a  horizontal  asymptote  as 
77      ->  00 .   Hence,  the  problem  can  be  transformed  to  an  initial  value 
problem  by  replacing  the  condition,  as  shown  in  equation  15,  by: 

f '  (0)  =  a>  *  0   .  [17] 

By  varying  the  parameter  00,  the  condition  of  equation  15  can  be  satisfied  for 
various  values  of  9Qi    provided  the  initial  value  problem  can  be  solved.   This 
solution  is  possible  by  employing  the  Taylor  series  method.   To  that  end,  the 
nth  derivative  of  equation  13  at  7?  =  7?  0  was  evaluated,  and  the  resulting 
equation  was  solved  for  the  derivative  of  order  (n  +  2): 

n 

f  (n+2)q)  -  -[n/^vv  *  ■*  (a)g0'  -  L  0  * (k)  gio)  t(n-yj  [18] 

f  (7Jo)        •  =  1,  2,  3,-  — 
The  Taylor  series  solution  is  of  the  form: 

« <H>  ■  &  « (k)Q|.>  (?M°)k 

k!  .  [19] 


The  coefficients  f 


(M 


ion  1Q.   Th 


were  computed  from  equation  l8,  and  f  (7?) 


was  evaluated  by  equation  19.   These  computations  were  made  on  the  IBM  65O. 


} 


Series  19  vill  have,  in  general,  a  finite  radius  of  convergence.   To 
extend  the  range  of  validity,  analytic  continuation  was  used  (k) .   In  this  way 
f  (7|)  was  delineated  for  7|  >  0.   Since  JJ (  f   ,  f   )  =  f  (f\) ,    equation  11 
becomes:  -* 


=  ^  In  f  (7|),  7\>   0 


[20] 


The  flux  J  at  the  boundary  x  =  0,  is  given  by: 


M 


X  =  0 


d  (0)  a© 

ax 


X  =   0 
By  the  use  of  equations  7,  17,  and  20  this  can  be  written: 


iW . 


[21] 


The  total  quantity  of  water  accumulated  inside  the  given  region  in 
time  t  is  given  by: 


Q  (t)  =  -  A 


which  in  view  of  equation  21  becomes: 


J 
,  J 


dt 


x  =  0 


Q  (t)  =  -  oA  V  2  0C  t 

jr  [22] 

Horizontal  infiltration  and  outflow  problems  are  characterized  by  co  <  0  and 
a)  >  0  respectively,  as  can  be  seen  from  equations  17  and  21. 

Linear  Diffusivity. 

The  analysis  in  the  case  when  the  diffusivity  function  is  linear  and 
of  the  form 

D  (0)  =  0(9  +  A  (  OC  >  0,  £  >   0)  [23] 

where  oc  ,  ^   are  parameters,  is  analogous  to  that  used  in  the  exponential  case. 


k/   The  use  of  this  or  other  patented  equipment  in  this  study  does  not  imply 

approval  of  the  product  to  the  exclusion  of  others  that  may  also  be  suitable, 


The  boundary  value  problems  to  be  solved  here  are: 

0  <  X  <00  ;  t  >  0, 

9  (0,       t)   =  9k,  (k  ::  1,   2),      t  >  0, 


d9  _  d 

"5t  "5c 


(oc  e+p  )    f 


[2U] 
[25] 


9  (x,  0)  =  9q,        x  >  0, 


0  <  9P  <  9  <  9, 


o    '1    .  [26] 

Again  k  =  1  characterizes  the  inflow  problem  and  k  =  2  the  outflow  problem. 
Solutions  of  equation  2U-  of  the  form 

y        M  J  ,  f  )  =  f  (7|)  [27] 

were  sought,  where  7?  =  ,.x=— -  ,  £  =  —  ,  and  ^*  =  ~"p  •   Equation  2U  becomes 

[f  (/J)  +A]  •  f"  (7J)  =  -  f  (7|)  [  7[  +  f  (7|)]  , 

where  A  =  ^/cC  .   Finally,  by  making  the  substitution 

f  (7J)  =  f  Oj)  +  A 

the  preceding  problems  take  the  form: 

F  (7])  F"  (7])  =  -  F"  (7J)   [  7\     +  F"  (7J)  ]  , 
F  (0)  =  9k  +A  , 
F  (©o)  =  9^  +A  , 


[28] 

[29] 

[30] 
[31] 


0  <  9P  <  9  <  9, 


[32 


By  evaluating  the  nth  derivative  of  equation  30  at  7?  =  7?   the 
following  relationship  can  be  obtained: 

n     ^  **  J 


+  .£  (")  /F<k+1» 


k=2<k>  <F 


(n-k+1)        (k)      (n-k+2)] 

a  =  2,  5,  h— 


[35] 


The  Taylor  series  solution  is  of  the  form: 

•  O  • 

(k) 
The  coefficients  F  )y%  \   were  computed  from  equation  35  and  F  (77) 


[5M 


was  evaluated  from  equation  3^-«   These  computations  were  made  on  the  IBM  65O. 

In  a  manner  similar  to  that  used  in  the  exponential  diffusivity  case, 
the  condition  of  equation  32  was  replaced  by: 


r-  (7?o)  = 


CO 


The    flux,    J,    at  x  =  0  and   the    total   quantity  of  water   accumulated  in 
the    region  are    given   by: 


tj]       -  - « (\ +  x  u  rsi 

X    =     0 


2t  [35] 


and 


Q   (t)   =  -  a)  (0k  +  A  )  A-\/2CXt        ^  6] 

Again,  to  <  0  denotes  the  inflow  problem  and  00  >  0  the  outflow 
problem. 

Method  of  Computation. 

The  procedures  for  obtaining  the  power  series  solutions  of  equations 
8  and  2k   for  the  exponential  and  linear  diffusivity  functions  were  similar. 
Flow  diagrams  outlining  the  computational  program  are  given  in  Appendix  III. 
The  program  consisted  of  the  following:   (a)  computation  of  the  Taylor  series 
coefficients  f  (k)  /n    ,  or  F  (k)  ,yi   \,   (b)  evaluation  of  the  Taylor  series  for 
f  (77)  or  F  (7?)  in  its  interval  of  convergence.   The  interval  of  convergence  was 
found  to  be  finite.   This  necessitated  extending  the  power  series  solution  to  a 
larger  interval  of  convergence  by  the  process  of  analytic  continuation  (h). 
When  oscillations  in  the  values  of  f  (7?)  or  F  (77)  were  detected,-^'  a  value  of 
7?=7L  just  to  the  left  of  the  right  endpoint  of  the  interval  of  convergence  was 
selected.   The  function  f  (77)  (or  F  (71))   and  a  sufficient  number  of  its  succes- 
sive derivatives  were  computed  at  '(1,  which  determined  a  Taylor  series  solu- 
tion defined  over  an  interval  of  convergence  77^  <  77  <  ^[2-      The  above  computa- 
tional procedure  was  continued  until  the  difference  between  successive  ordinates 
was  less  than  some  preassigned  quantity. 

A  particular  pair  of  values  of  0^  and  a>  were  selected  and  the  compu- 
tational procedure  described  above  was  carried  out.   Corresponding  to  the  0^ 
and  cd  chosen,  a  particular  value  of  9Q   was  obtained.   By  varying  cu  with  a  fixed 


_£/  A  series  consisting  of  2h   terms  was  used  to  evaluate  f  (77)  or  F  (77) 

8 


Gu   a  series  of  problems  was  solved.   These  computed  results  are  shown  in 
figures  1,  2,  3,  and  k.      The  date  in  all  the  other  figures  in  this  report  are 
based  on  these  computed  results  as  explained  below. 

Results  -  Exponential  Diffusivity. 

The  water  content,  0,  is  related  to  f  (71)  by  the  equivalent  equations 

f  (J|)  =  e  /39   or  In  f  (7/)  =  fi&        .  [37] 

In  order  to  use  the  computed  results  in  practice,  transformations  must  be  made 
into  measurable  quantities.   For  that  purpose,  the  following  water  content 
ratios  are  convenient  to  use: 

e  -  qp  -  m  f  (7?)  -  m  f  (o) 

9o  -  92  "  in  f  M  -  In  f  (0)  '  °UtUOW  [58] 

9  -  So  _  In  f  (7?)  -  In  f  (oo) 

Q  -   G     "  In  f  (0)  ■  In  £  W  [39] 


Oa  6  a. 

where  f  (0)  =  e   2  in  the  outflow  case,  f  (0)  =  e    *  in  the  inflow  case,  and 
f  (oo)  =  e  P^o   in  both  cases.   In  this  way  these  ratios  will  have  values  in  the 
range  from  0  to  1.0.   From  the  results  of  the  computations  of  f  (7?)  for  various 
values  of  the  parameter  oo,  In  f  (7?)  was  computed  for  these  same  values  of  the 
parameters  and  the  ratios  in  J>Q   and  39  were  obtained.   It  was  more  convenient 


to  graph  (G  -    G2)/(G     ■  92)  as  a  function  of  the  ratio 


x  /v^r=7?  V2/  e^ 


[^o; 


and   (Q   -   G    )/(G.    -    Q    )    as   a   function  of   the   ratio 
o  1  o 


x    /-yrTt"  =  7?  -^2/  e  -  -o 


/3  Q 

[kl] 


Figures  5,  6,  and  7  show  the  plots  for  inflow  and  outflow,  respectively,  for 
various  values  of  the  parameter  cu  =  f  (0). 

The  procedure  for  translating  the  dimensionless  quantities  of  figures 
5,  6,  and  7  into  the  corresponding  specific  dimensional  quantities  of  a  given 
problem  is  as  follows:   For  a  given  porous  medium  the  diffusivity  function  must 
be  known.   The  initial  water  content,  QQ,    and  the  water  content  at  x  =  0,  0^ 
must  also  be  given.   When  k  =  1,  i.e.  infiltration,  the  ratio  D^/D0  is  equal  to 

eP\l  o)  where   0j_  >  G0.      When  k  =  2,    i.e.    outflow,    the   ratio   D0/D2   is   given 

by  e  °  2)  where   QQ  >  62 .      Since    the   diffusivity   function   is   assumed   to  be 

exponential    and  known,    both  oc  and /3  are   known. 
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FIGURE  1. --Original  computed  results  for  inflow  with  exponential  diffusivity. 
The  parameter  on  the  curves  is  oj  =  f  (0). 
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EXPONENTIAL  -  OUTFLOW 


IOOO 


FIGURE  2. --Original  computed  results  for  outflow  with  exponential  diffusivity. 
The  parameter  on  the  curves  is  go  =  f '  (0). 
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LINEAR    DIFFUSIVITY- INFLOW 


100 


F(iy)      - 


FIGURE  3. --Original  computed  results  for  inflow  with  a  linear  dif fusivity, 
The  parameter  on  the  curves  is  co  =  F1  (0). 
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LINEAR    DIFFUSIVITY- OUTFLOW 
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FIGURE   5.--Dimensionless  plot  of  water  content  versus-|fD0t   for   inflow  with  an 
exponential    diffusivity.      The   parameter   is  u)*. 
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To   determine    the   value   of  co  corresponding   to   given  values   of  0O  ^  0 
and  9i    the    following  procedure    is   used  in  connection  with    the   exponential 
diffusivity  inflow  case.      Compute    the   ratio   D^/Dq    from  e  P^l/e    "®o.      Find  a 
value   of  co*  corresponding   to   this   diffusivity  ratio  by  use   of    figure  8. 
Calculate  co  from: 


CO  =  CO* 


This  equation  is  derived  in  Appendix  I. 
figure  8  is  discussed  in  Appendix  II. 


'\fJT°    • 

The  construction  of  the  graph  in 


[U2] 


The  outflow  case  is  treated  in  a  similar  manner.   The  ratio  Do/D2  is 
computed  from  e  ^  9o/e  ^  ®2   and  a  value  of  co*  found  by  the  use  of  figure  8. 
The  appropriate  value  of  co  corresponding  to  the  given  0g  is  then  obtained  from: 


CO  =  CO* 


V- 


0    9, 


[*3] 


The  flux  and  accumulated  water  content  can  now  be  found  from  equations  21  and 
22. 

The  magnitude  of  the  ratio  of  the  outflow  flux  to  the  inflow  flux 
when  (a)   l#k  -  9Ql       is  the  same  for  outflow  as  for  inflow,   (b)  when  92    for 
outflow  is  equal  to  9Q   for  inflow,  and   (c)  when  0^  for  inflow  is  equal  to  9Q 
for  outflow  is  given  by:  , — 

a)  outn  7°^ 
V  2t 


[J]  o"t 


[J] 


in 


x  =  0 


__£ 


co  in 


co  out 

CO*    outj 

co  in 

co*   in  [ 

2t 


Under  the  above  conditions  D0/D2  =  D^/Dq.   From  figure  8  and  by  assuming  the 
above  conditions,  note  that  the  inflow  flux  is  larger  in  magnitude  than  the 
outflow  flux. 

The  results  of  the  computations,  as  shown  in  figures  5*  6,  and  7> 
could  be  presented  with  a  diffusivity  ratio  as  the  parameter,  as  Gardner  (l) 
has  done.   However,  since  the  computations  were  originally  made  with  co  and  9\^ 
as  the  boundary  conditions  it  was  convenient  to  use  co*  as  the  parameter  in 
the  curves  of  figures  5,  6,  and  7«   Also,  if  the  flux  is  desired,  the  quantity 
co  is  needed. 

Results  -  Linear  Diffusivity. 

The  water  content  is  related  to  F  (J))   by 


(7|)=e  +  A, 


[kk] 

where  9  is  defined  in  27  and  A  =  ft  I  &-  .  The  graphs  for  the  inflow  and  outflow 
cases  are  shown  in  figures  9  and  10,  respectively,  for  the  specific  case  /\  =  0. 
It  was  convenient  to  plot  the  ratio  9/9-^   versus  x/  ~]j   2D-,  t  in  the  inflow  case 
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FIGURE  8. --Graph  of  diffusivity  ratio  versus  ay*.   The  construction  of  this 
graph  is  explained  in  Appendix  II. 
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FIGURE   10.--Dimensionless  plot  of  water  content  versus  x/  y2D2t   for   outflow 

with  a  linear   diffusivity   function  in  which    ^    =  0.      The   parameter 

is    03*. 
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and  the  ratio  0/02  versus  x/"V  2D2t  in  the  outflow  case.   In  each  graph  the 
curves  are  labeled  with  co  as  a  parameter. 

The  procedure  for  obtaining  the  appropriate  value  of  co  for  a  given 
pair  of  values  of  0O  and  0^  which  define  a  specific  flow  problem  is  as  follows. 
Calculate  DQ  and  D^  by  the  use  of  equation  2J>.      Form  the  ratio  D0/D2  of  D^/D0 
and  use  figure  11  to  obtain  or*.   This  assumes  X     -    0.   Calculate  co  from 


cd  =  co*"V  0^ 

Next  select  the  appropriate  curve  in  figures  9  or  10  which  is  labeled  with  the 
above  value  of  co. 

The  flux  and  accumulated  quantity  of  water  can  be  calculated  from 
equations  35  and  36. 


by: 


The  magnitude  of  the  ratio  of  inflow  flux  to  outflow  flux  is  given 


[J]  in 


[J]  out 


x  =  0 


oo  in   (01  +  A  ) 


CO 


OUt  (02  +  )  ) 


[^5] 


If  we  assume  (a)  that  0^  for  inflow  is  equal  to  0Q  for  outflow  and  (b)  that  02 
for  outflow  is  equal  to  0O  for  inflow,  then  from  the  equation  for  the  diffusiv- 
ity,  equation  2J,  it  follows  that 


'2I| 


IN 


JDft 
D2 


(D0)  IN  t   (Dj  OUT 


OUT 


Inspection  of  figure  11  shows  that  for  a  given  diffusivity  ratio  [co  in|  is  less 
than|co  outj. 


If  A  =  0  then: 
J  in 


J  out 


x  =  0 


to  in  01 


CO 


OUt  02 


">*  in_   /0i^ 


<B*    OU 

co*  in 


out 


5/2 1 

DM  3/2 
D°  in 


iDa\3/2 


to*  out 


'2  out 


The  magnitude  of  this  ratio  can  be  computed  with  the  aid  of  figure 
11.   The  computations  show  that  the  inflow  flux  is  larger  in  magnitude  than  the 
outflow  flux.   Some  computed  results  are  shown  in  Table  1. 
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FIGURE  11. --Graph  of  diffusivity  ratio  versus  oj*.   The  construction  of  this 
graph  is  explained  in  Appendix  II. 
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TABLE 

1. 

Dl 
Do 

in 

or 

Do 

D2   out 

J  in 
J  out 

k 

1.12 

8 

1.24 

12 

1.29 

Applications. 

If  a  flow  system  can  be  described  by  equation  1  and  if  the  boundary 
conditions  stated  in  equations  9  and  10  are  applicable,  then  it  is  possible  to 
predict  the  water  content  as  a  function  of  time  and  position  in  the  flow 
system.   The  flux  and  the  total  amount  of  water  accumulated  can  also  be  calcu- 
lated.  If  the  diffusivity  is  assumed  to  be  exponential  or  linear  the  constants 
in  the  diffusivity  function  can  be  obtained  from  a  measurement  of  the  water 
content  distribution  at  a  given  time,  the  appropriate  computed  distribution  of 
water  content  as  shown  in  figures  5,  6,  7*  9,    or  10,  and  the  value  of  the 
diffusivity  at  one  water  content. 

The  diffusivity  data  for  two  soils,  Sarpy  loam  and  Geary  silt  loam, 
appear  to  be  of  exponential  nature.   We  have  determined  the  water  content  as  a 
function  of  distance  for  several  different  boundary  conditions  as  shown  in 
figures  12  and  13.   For  Sarpy  loam  for  infiltration  with  9Q   =  0.06l  and 

.9l  =  0.U1  the  ratio  Dl/D0  Is:   DX/D0  =  e  *<81  "  9°>  =  e  ^  «>-M   =  *»■ 

From  figure  8,  for  inflow,  03*  is  seen  to  be  about  -79«5«   By  using 
equation  k2   00  is  calculated  to  be  -Ik1?.      The  distributions  of  water  content 
shown  in  figure  12  were  calculated  from  the  curves  of  figure  5  labeled  with 
oj*  =  -79«67«   ^he  other  curves  in  figure  12  were  obtained  in  a  similar  fashion. 

The  flux  at  x  =  0  for  Sarpy  loam  with  0^  =  O.Ul  and  9Q   =  0.06l  is 
given  by  equation  21:  , r- 

r      n  -Ikk         "A   I  l.2k    X     10 


x  =  0  I9.78       V  2t 

s  3-73  X  10'2 

-yTt  (cm/ sec) 

The    total   quantity  of  water  accumulated  in  time    t   is   given  by  equation  22: 

Q    (t)   =  -  j^£       V2  X  1.24  X  10"4  X  t 

=  0.115  k~\ft  (cm5) 

Figure  13  shows  9   versus  x  for  indicated  values  of  t  for  outflow. 
These  curves  were  constructed  in  a  manner  similar  to  that  described  above. 

Note  that,  for  the  exponential  diffusivity  case,  the  water  content 
scale  can  be  translated  up  or  down  without  influencing  the  shape  of  the  (0vs  x) 
curve.   That  is,  the  data  of  figure  13  could  apply  to  any  boundary  conditions 
where  (9\    -    9Q)    is  constant.   For  example,  the  data  for  Sarpy  loam  with 
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SARPY  LOAM     a  =  1.24  X I0'4  cm2/sec 
/?=  19.78 
D=ae#* 


253  sec 


j_ 


J_^ L 


&0=<273 


0O  =0.061 


GEARY  SILT  LOAM 


—   0,= 0.386 


a  =  5.66  XIO"6 
fi=  21.0 
D=  ae^0 


6020  sec 
1720  sec 
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0O=  0.053 
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FIGURE  12. --Water  content  versus  distance  for  horizontal  inflow  for  Geary  silt 
loam  and  Sarpy  loam  vith  boundary  conditions  as  indicated. 
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■Water  content  versus  distance  for  horizontal  outflow  for  Geary  silt 
loam  and  Sarpy  loam  with  boundary  conditions  as  indicated. 
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9\   =  O.^l  and  9Q  =   0.06l  would  also  apply  to  the  case  with  9\   =   0.460  and 

9Q   =  0.111.   The  ratio  D>1  would  be  the  same  but,  of  course,  the  flux  would  be 

Do 
changed,  because  go  would  be  changed.   The  value  of  a.  corresponding  to 

01  =  0.46  and  0O   =  0.111  is 

„^v/   19.78  *  0.111 

=  -79  X  3.O3 

=  -239 

This  value  of  a.  would  be  used  to  compute  the  flux  for  these  new  boundary 
conditions. 

Similar  application  could  be  made  of  the  linear  data  if  the  linear 
diffusivity  function  is  known  and  if  the  boundary  conditions  agree  with  those 
of  equations  25  and  26. 

Comparison  of  Linear  and  Exponential  Cases. 

In  order  to  make  a  comparison  between  an  outflow  system  in  which  a 
linear  diffusivity  function  applies  and  one  in  which  an  exponential  function 
applies,  we  will  impose  the  conditions  that  (a)  for  9  =  92    the  exponential  and 
linear  dif fusivities  are  the  same  and  (b)  for  9   =  9Q   the  exponential  and  linear 
dif fusivities  are  equal.   These  conditions  are  shown  graphically  in  figure  Ik. 

x 

The  water  content  as  a  function  of  ~^/2D2t   is  shown  in  figure  15  for 
two  examples  where  92   has  been  set  equal  to  0.01  and  90   equal  to  0.106  and 
0.21.   It  is  evident  from  this  figure  that  the  exponential  case  does  not  drain 
as  rapidly  as  the  linear  case  under  the  imposed  conditions. 

The  outflow  flux  for  the  linear  case  is  given  by  application  of 
equation  35  to  De: 

J  =  o.318-v/d^~ 

where  the  appropriate  to  and  OC values  have  been  determined  from  the  diffusivity 
ratio  and  92.      The  outflow  flux  for  the  exponential  case  is  given  by  applying 
equation  21 : 

J  =  0.223-y/l^ 

where  to,  j3   andOChave  been  determined  from  the  diffusivity  ratio  and  92«   From 
these  two  equations  it  is  seen  that  at  any  given  t  the  flux  for  the  linear  case 
is  larger  than  that  for  the  exponential  case. 
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FIGURE  lU. --Graph  of  D  versus  9   to  show  how  the  comparison  of  exponential  and 
linear  outflow  cases  shown  in  figure  15  was  made. 
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FIGURE   15, 


■Comparison  of   linear   and  exponential    outflow  cases, 


28. 


Estimation  of  Dif fusivity  Functions. 

If  it  is  assumed  that  the  diffusivity  function  is  of  either  the 
exponential  or  linear  type  and  if  the  value  of  the  diffusivity  at  the  water 
content  9Q   is  known,  the  constant  in  the  diffusivity  function  can  be  determined 
from  a  measurement  of  the  water  content  distribution  at  a  given  time  in  either 
an  inflow  or  outflow  system.   By  assuming  for  the  moment  that  the  diffusivity 
function  is  exponential,  that  DQ  is  known  and  that  measurements  are  made  on  an 
inflow  system,  the  resulting  9   versus  x  data  could  be  plotted  in  the  reduced 
form  of  figure  1.   By  matching  curves  and  interpolation  a  value  of  a>   could  be 
estimated  and  from  this  and  0O,  0]_,  and  D0  a  value  of  D^  could  be  obtained. 
Since  D0  and  D]_  are  now  known,  0  andod  can  be  determined.   A  similar  procedure 
could  be  followed  if  a  linear  diffusivity  is  assumed. 
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APPENDIX  I 


Derivation  of  Slope  Conversion  Formulas  for  the  Exponential 
and  Linear  Diffusivity  Functions 

(A)  Inflow  Case  for  the  Exponential  Function  f  (7|)  =  e     ( /&  >  0) .        =   X/l/Jctt 

The  plot  of  (0  -  0o)/(0i    -  9Q)   versus  x  /~[fDQt   gives  a  single  curve 
for  all  boundary  conditions,  where  (9\    -  9Q)  =   const,  or  D^/D0  =  const.   Also, 
7[  =  £  -\J   Do/(2  0C   ),  where  A  =  x  /-j/D0t. 


Now, 


co  = 


df 


(Tjl 


n 


)=  OL 


df    (Tp      d0   d4 
d9  dr   dTj 


&  e  d0 


=   0 


[yD°/(2  "7  d^ 


0O' 


(=0 


=V2  /Se  '(«i  -f)    141 


since   D     =  o£  e        °. 


*.. 


Similarly,  for  f*  (7J)  =  e  *  0*, 


CO* 


since  D0*  =  OC  e 


£  0rt* 


2  ; 


d9* 

4.0 


Moreover,  for  two  different  boundary  conditions  such  that 

the  two  curves  corresponding  to  co  and  co*  plot  on  one  curve  in  the  graph  of 
(d   -   0O)/(01  ~  0O)  versus  {    .   Therefore, 


[1] 


de 


d9*i 

d/ 


l%o    "l%o      . 


[2] 


Forming   the   ratio   of  co*   to  co  and  making  use   of  conditions   as    shown  in  equations 
1   and  2,    the    following  result   obtains 


If  9Q*  =   0  as   a   standard,    then  equation  3   becomes 

co  =  co*  Ve   *  9°        . 


[3J 
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(B)   Outflow  Case    for   the  Exponential   Function. 

Just  as   in   the   inflow  case,    the    following  expressions   are    found: 

1%  J/"-  o 


CO 


-T/i><.    KbllfcZ)      ]$ 


c 


By  making  use  of  the  relations 

D2*  =    OL  e 


&    92*,      D2   =   oCe 


0es 


'd9" 


de 


results   in   the    formula 


tA 


a>*  =  a)     \  /  e 
For  92*  =  0,    equation  5  becomes 


a)  =  cx>*\e 


/*(92*   -    92) 


£92 


[5] 


[6] 


(C)    Inflow  Case   for    the   Linear   Diffusivity  Function  D  =    PC  9  +  fij  (  CX-    >  0, 
/9  >  0),   7}  =  x    /y2ac  t. 

By   definition,    F   (7?)   =  9  +  X    ,   where  A  =    ^  /CC    and   D1    =0-9!   +^    . 
Also,    F   (0)   =  9L  +  X  and  7)  =  [  y5j7ot    ,   where  /  =  x   /-[/2Dlt. 


Now, 


0)    = 


dF 


ial 


fdF 


■  o-[ 


di 

£=< 

F   (.0)    I'd   (F    fW)| 
F    (0) 


r^  xk-  o 


V^£"   [<^L  o    "l/ntt  L^i-  -'«*=  o 


^fe    (ft)] 


f    IMO)^ 


=   0 


In   the   same  way, 


CD*   =    "l/F*    (0)   Ij- 

sincefL  urn        _fd_  /  f*  ffl) 

Since  dr     F  (o;     ,    _  n  -  (d/-     (F*  (0) 


F*    (7?) 

F*   (0 


£  ,0      [%     V*  (0)/J^o 


<=° 
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it   follows    that 


co    "  Vf     (0)     "Vol     +  A 


A  '  [7] 

In  the  report,  the  case  j3  =   0  or  J^=   0  was  treated.   Hence  the 
conversion  formula  for  this  case  is 


cd*  =  arjQi*   01    .  [8] 

In  particular,    if  Q\*  =  1, 

cu  =  CD*Ve7       .  [9] 

(D)  Outflow  Case  for  the  Linear  Diffusivity  Function. 

Here,    D2   =  oC92  +  0   ,   7}  =  £  ~]f~D^T<K   f   where    £  =  x/Y2D2t;    and 
F   (0)   =  92  +A    •      As   before  . 

n  /f*  (o)       n&*+A 
^  =  a)VF(o)    =aiVe^    +A       .  [L0] 

For  the  case  treated  here,  A  =  0,  so  that 


CD*  =  CD  l/92*  /  02    .  [11] 

In  the  particular  case  02*  =  1, 

cd  =  a>*yij~"  .  [12] 
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APPENDIX  II 

Construction  of  Graphs  in  Figs.  8  and  II 

(A)  Inflow  Case  for  the  Exponential  Function. 

Using  the  original  data,  corresponding  to  a  given  boundary  water 
content  9   =  0^**  ^  0  and  a  given  slope  a)  =  ay**,  an  initial  water  content 
0  =  0O**  was  determined.   The  values  9  =  9\**   and  9   =  9Q**   thus  determined  the 
ratio 

Dl     ,  fi   (01**  -  0O**)  _  f  (0) 

Do  =  f  W    •  [1] 

By  varying  the  slope  cjd  =  ox**,  various  ratios  D^/D0  were  computed.   However,  in 
view  of  the  fact  that  the  ratio  Dl/D0  is  the  same  whenever  the  difference 
©1**  -  0O**  =  constant  (e.g.,  for  0^**  =  0.20,  9G**  =  0.10  as  well  as  for 
01**  =  0.10,  0O*"*  =  0.00),  it  was  decided  to  use  zero  boundary  water  content  in 
this  case  as  a  standard  reference  point.   To  determine  the  corresponding  change 
in  the  slope  ay,  appeal  was  made  to  formula  3  of  Appendix  I. 

If  ay**  corresponds  to  0O**  and  o>*  corresponds  to  0O*,  this  formula 
reads  

a,*  =  a**-\/e    ^"°*-    9°">        .  [2] 

Therefore,    for   0O*  =  0,    formula  2  becomes 


hJP^' 


03*  =  ay**     /\ler~  wo**        .  [3] 


It  was  this  formula  which  was  used  to  draw  the  graph  in  Fig.  8  of  ay*  versus 
Di/Do. 

Suppose  a  slope  a>  is  desired  corresponding  to  any  given  values 
9  =  9\   and  9   =  9Q.      Then,  these  values  determine  the  ratio  Di/D0  which  in  turn 
determines  ay*  corresponding  to  0O*  =   0  from  Fig.  8.   The  new  adjusted  slope 
now  follows  from  formula  k,   Appendix  I,  namely, 

u)  =  03*lje  P  Qo        .  [k] 

(B)  Outflow  Case  for  the  Exponential  Function. 

Using  the  original  data,  corresponding  to  a  given  boundary  water  con- 
tent 9   =  02**  ^  0  anc*  a  given  slope  ay  =  ay**,  an  initial  water  content  9  =   0O** 
was  determined.   The  values  9   =  02**  and  9   =  0O**  thus  determined  the  ratio 

Do      P  (90**   -  02**)    j    (°*) 

D2  "  e  ~  f  (0)    .  [1] 


3h 


By  varying  the  slope  co  =  co**,  various  ratios  D0/D2  were  computed.   However, 
in  view  of  the  fact  that  the  ratio  D0/D2  is  the  same  whenever  the  9Q**  -   92** 
=   constant,  (e.g.,  for  9Q   =  0.20,  9\    =  0.10  as  well  as  for  9Q   =  0.10, 
92   =  0.00),-  it  was  decided  to  use  zero  boundary  water  content  in  this  case  as 
a  standard  reference  point.   To  determine  the  corresponding  change  in  the  slope 
co,  appeal  was  made  to  formula  5  of  Appendix  I.   If  co**  corresponds  to  92**   and 
to*  corresponds  to  92* ,  this  formula  reads 

e  '  v  2     2   ;  [2] 


Therefore,  for  02*  =  0>  formula  2  becomes 


/3  92 


** 


CO*  =  to**  I\je   ^    Z  •  [3] 

It  was  this  formula  which  was  used  to  draw  the  graph  in  Fig.  8  of  to*  versus 
D0/D2- 

Suppose  a  slope  co  is  desired  corresponding  to  any  given  values 
9  =   92   and  9  =  9Q.      Then,  these  values  determine  the  ratio  D0/pp.  which  in  turn 
determines  co*  corresponding  to  92*  =   0  from  Fig.  8.   The  new  adjusted  slope  now 
follows  from  formula  6,  Appendix  I,  namely, 


4' 


co  =  co*\/e   "  U2        .  [1+] 

(C)  Inflow  Case  for  the  Linear  Diffusivity  Function  D  =  OC  9    (  /3   =  Q). 

From  the  original  data,  corresponding  to  a  given  boundary  water 
content  9   =  9\**  £   0  and  a  given  slope  co  =  co**,  an  initial  water  content 
9   =  9Q**  was  determined.   The  values  9   =  9i**  and  9   =  9Q**   thus  determined  the 
ratio 

RL  9l**     F  (0) 

Do  ~  0O**  ~  F  H    .  [1] 

By  varying  the  slope  co  =  co**,  various  ratios  Di/D0  were  computed.   Here,  in 
view  of  equation  1,  a  given  ratio  Di/D0  =  constant  will  correspond  to  all 
problems  in  which  the  ratio  9\**/90**   is  equal  to  this  same  constant.   In  this 
case,  since  a  zero  water  constant  as  a  base  could  not  be  used,  it  was  decided 
to  base  co  on  a  water  content  of  unity. 

From  formula  8,  Appendix  I, 

co*  =  co**  V  01*/9i**   .  [2] 

For  9\*   =  1,  this  formula  becomes 

co*  =  co**/^**    .  [3] 

The  graph  in  Fig.  11  was  constructed  on  the  basis  of  this  formula. 
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Suppose  a  slope  co  is  desired  corresponding  to  any  given  values 
9   =  9\   and  9  =   0o.   Then,  these  values  determine  the  ratio  D^/D0  which  in 
turn  determines  co*  corresponding  to  01*  =  1  from  fig.  11.   The  new  adjusted 
slope  co  now  follows  from  formula  9*  Appendix  I,  namely, 

co  =  (x>*"\J9i 

(D)  Outflow  Case  for  the  Linear  Diffusivity  Function  D  =  (X  9   (  fi    =  0). 

In  this  case  the  ratio 

Do.    90**  F  (oo) 

D2  ~  e2**  ~  F  (0) 

was  computed.   As  for  inflow,  unity  was  taken  as  the  base  for  computing  to. 
From  formula  11,  Appendix  I, 


co*  =  co**  Ye2*/e2**   , 

which  for  92*   =  1  becomes 


7  Ye 


to*  =  to**/^e2** 

Formula  3  was  used  to  construct  Fig.  11. 

The  slope  co  corresponding  to  any  given  values  9  =  92,    9   =  9o   is 
obtained  as  follows.   Using  the  preceding  values,  determine  the  ratio  Dq/E^. 
From  fig.  11  obtain  the  value  of  co*  corresponding  to  02*  =  1.   The  new  adjusted 
slope  co  now  follows  from  formula  12,  Appendix  I, 


to 


=  co*!^" 
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APPENDIX  III 


Flow  Diagrams  Outlining  the  Computational  Program 
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Enter  Here 


■F|  =  Pitt  A,  F2=  Pitt  B 
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EN 
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no 


no 


no 
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LINEAR     DIFFUSIVITY    SUBROUTINE 


n-k  +  2=2,  EN  =  I,  n  =  2 
F3  =  F2(F2+77o);    F4  =  (F34I)F2+(2F2+t7o)F3 
-F,                                       -F, 

1 

! 

k  =  2,EN=EN  +  l  ,KP=2,A=0,  NM  =  EN-EN',   NN  =  NM/2 

k'  =  k  +  l 
NN'  =  NN,NM 

KP' 
KP'=  KP  +  I 

NM'  =  NM-I 

<! 

no 

A'=    NN  (Fk+lFn-k  +  ,+Fk  Fn"k+2)+A 

n'  =  n  +  l 

?  EN  =  KP  ? 

il 

yes 

F 

n+2=[(EN  +  l)F2+i70]Fn+l  +  (l  +  F3)nFn  +  A 

I 

-F| 

?     n=2l  ? 

yes 

-® 

' 

r" 

EXPONENTIAL    DIFFUSIVITY    SUBROUTINE 

T 


n  +  2-k=2,n  =  1 

1 

' 

EN  =  EN+I  =  NN  =  NM,  KP=l,k  =  l,  A  =  0 

k'=k  +  l 
NN'=NN,NM 
KP' 

KP'=KP  +  I 

NM'=NM-I 

! 

1 

no 

NNFkFn-k+2+A  =  A- 

EN  =  KP 

i 

,yes 

Fn< 

h2  =  nFn  +  i70Fn  +  l  +/ 

* 

"F, 

n=22 

no  __ 

n'=n  +  l 

yes 

D 
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